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Abstract. We present the constraints on the cosmological parameters obtained with the Ep^^ic-E^ correlation found with the 
most recent sample of 19 GRBs with spectroscopically measured redshift and well determined prompt emission spectral and 
afterglow parameters. We compare our results obtained in the two possible uniform jet scenarios, i.e. assuming a homogeneous 
density profile (HM) or a wind density profile (WM) for the circumburst medium. Better constraints on Q.m and SIa are obtained 
with the (tighter) i5pcak-i57 correlation derived in the wind density scenario. We explore the improvements to the constraints 
of the cosmological parameters that could be reached with a large sample, ~ 150 GRBs, in the future. We study the possibility 
to calibrate the slope of these correlations. Our optimization analysis suggests that ~ 12 GRBs with redshift z £ (0.9, 1.1) 
can be used to calibrate the -Epcak-S^ with a precision better than 1%. The same precision is expected for the same number of 
bursts with z G (0.45, 0.75). This result suggests that we do not necessarily need a large sample of low z GRBs for calibrating 
the slope of these correlations. 
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1. Introduction 

Gamma Ray Bursts (GRBs) are presently detected out to very 
high redshifts (the new limit being GRB 050904 at z = 6.29, 
^ Kawai et al. 2005) and this makes them extremely attractive for 
■ observational cosmology. They might have profound impact on 
' (i) the study of the epoch of reionization, (ii) the characteriza- 
• tion of the properties of the cosmic intergalactic medium, (iii) 
' the description of the geometry of the Universe and (vi) the 
investigation of the nature of Dark Energy. 

GRBs have recently been discussed (Ghirlanda et al. 2004b, 
Firmani et al. 2005, Lamb et al. 2005) as a possible new class of 
"standard candles" to be used to constrain the Universe models 
(Ghirlanda et al. 2004b, Firmani et al. 2005, Liang & Zhang 
2005). What "standardize" the GRB energetics is a tight rela- 
tion between their rest frame collimation corrected energy E-^ 
and the peak energy i?pcak of their vF„ prompt emission spec- 
tra (Ghirlanda et al. 2004). 

In the "standard" GRB scenario the collimation corrected 
energy is defined as — E^^\so{^ — cos6'j), where S-y.iso is 
the isotropic equivalent energy and 9-^ is the jet opening angle. 
This parameter can be derived from ij, i.e. the time when the af- 
terglow light curve steepens, in two different scenarios: (a) as- 
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suming that the circumburst medium is homogeneous (HM) or 
(b) assuming a stratified density profile produced, for instance, 
by the wind of the GRB progenitor (WM). In both cases the jet 
is assumed to be uniform. 

Originally the E-p^ak-E^ correlation was derived in the HM 
case with 15 GRBs (Ghirlanda et al. 2004) and resulted in a 
small scattered correlation (i?pcak £^°-^ with a ^ 0.15) 
which was used to derive interesting (though shallow) con- 
straints on the cosmological parameters (ilM, ^^a) and on the 
equation of state parameters of the Dark Energy (Firmani et al. 
2005). 

Recently Nava et al. 2005 (N05) derived the Epc!i\i-E^ cor- 
relation in the WM case. This correlation is less scattered, i.e. 
a = 0.08, than in the HM case and it is linear, i.e. -Epcak oc E^. 
N05 discussed its implications for the understanding of the dy- 
namics and radiative processes of GRBs. Although the i?pcak- 
E^ correlations (both in the HM and WM case) are model de- 
pendent, their consistency with the completely empirical cor- 
relation between ii^poak, and found by Liang & Zhang 
(2005 - hereafter LZ05), suggests that the model parameters 
are not much dispersed. It is therefore worth to derive the cos- 
mological constraints with the Ghirlanda correlation in both the 
HM and WM scenario and compare the results. 

In order to use GRBs as cosmological tools, through the 
above correlations, three foundamental parameters, -Epcak, E^ 
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and tjct, should be accurately measured or inferred. This re- 
quirement also applies to the empirical correlation of LZ05. 
Therefore only a limited number of GRBs, i.e. 19 out of ^ 
70 (up to Nov. 2005) with measured z, can be used as stan- 
dard candles. It also appears that the limited energy range of 
BAT onboard Swift (15-150 keV) allows to constrain only with 
moderate accuracy the i?pcak of particularly bright-soft bursts. 
However, given the perspective of the cosmological investiga- 
tion through GRBs, it is worth exploring the power of using 
GRBs as cosmological probes. 

Another still open issue related to the use of GRBs as 
standard candles is the so called "circularity problem" (see 
Ghisellini et al. 2005). This is due to the fact that the small 
number of GRBs with spectroscopic measured redshifts and 
their wide dispersion in z does not allow to calibrate the 
Ghirlanda correlation, which, in turn, depends on the cosmo- 
logical parameters that we want to constrain. To the aim of cal- 
ibrating this spectral-energy correlation, GRBs at low redshift 
are required. In fact for 2; < 0.1 the difference in the luminos- 
ity distance computed for different choices of the cosmologi- 
cal models [for JImj^^a G (0, 1)] is less than 8%. However, 
if (long) GRBs are produced by the death of massive stars, 
they should roughly follow the cosmic star formation history 
(SFR) and we should expect that the rate of low redshift events 
is quite small at z < 0.1. Instead, a considerably large number 
of GRBs with z > 1 should be collected in the next years by 
presently flying instruments (Swift and Hete-II). At such large 
redshifts, the cosmological models starts to play an important 
role. However, if it will be possible to have a sufficient number 
of GRBs with a similar redshift, it might still be possible to cal- 
ibrate the slope of these correlations with high redshift GRBs. 

The aim of this work is to find the cosmological constraints 
by using the Speak — £'7 correlation in the HM and WM case 
(Sec. 3) and to describe the level of precision that could be 
achieved with a large population of bursts (Sec. 4). Finally, we 
study the possibility to calibrate the slope of these relations 
with GRB samples at intermediate redshifts (Sec. 5). 

2. The correlation in tlie HM and WM 

case 

The jet opening angle 6'j in the two HM or WM density profile 
scenarios is: 

/ \ 3/8 / ^ ^ N 1/8 
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where n is the constant circumburst density in the HM case 
(Sari et al. 1999) and A* is the normalization of the density pro- 
file n(r) = 5 X lO^^A^r^^ g cm^'^ in the WM case (Chevalier 
& Li 2000). 

N05 derived the E^^^^^-E^ correlation in the HM and WM 
case with 18 GRBs (sample updated to Sept. 2005) for which 
secure measurements of z, i?pcak and have been reported 
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Fig. 1. Constraints on the cosmological parameters f^M, i^A 
obtained with the updated sample of 19 GRBs presented in 
Tab. 1 of N05 (to which GRB 051022 was added) in the ho- 
mogeneous density case (HM). The solid (red) contours, ob- 
tained with the 19 GRBs alone, represent the 68.3%, 90% 
and 99% confidence regions. The center of these contours (red 
cross) corresponds to a minimum = 15.25/17 dof and has 
= 0.23 and JIa = 0.81. The contours obtained with 
the 156 SN la of the "Gold" sample of Riess et al. 2004 are 
shown by the dashed (blue) lines. The joint GRBh-SN con- 
straints are represented by the shaded contours. We also show 
the 90% confidence contours obtained with the WMAP data 
(from Spergel et al. 2003). 

in the literature. While completing this work the Hete-II and 
Konus-Wind sateUites detected a new GRB (051022 - Olive et 
al. 2005) for which all the observables (i.e. z — 0.8 [Butler et 
al. 2005], £^pcak, ^j) required to compute the E^^^y-E-f correla- 
tion were measured. We therefore added this burst to the sam- 
ple of 18 events of N05. The spectrum of GRB 051022 as ob- 
served in the widest energy range by Konus-Wind (Golenetskii 
et al. 2005) is fitted by a powerlaw model with an high energy 
exponential cutoff: the low energy photon spectral index is a = 
-1.17 ± 0.04 and the {EFe) peak energy is £;pcak=510±35 
keV. The 2keV-20MeV energy fluence is 2.61±0.08xlG-4 
erg/cm^. The analysis of the X-ray light curve has shown the 
presence of a possible jet break at tb = 2.9±0.2 days (Racusin 
et al. 2005). For a standard cosmology with il^/ = 0.3 and 
= h = 0.7 (and assuming typical values for the den- 
sity and the wind profile normalization - see N05), we derive 
6'j(HM)=6.37±0.75° and 6'j(WM)=3.3 ± 0.1° for the HM and 
WM case. With the addition of this new bursts to the sample of 
18 events of N05, we find Spoak oc £:0-67±o.o4 (^j^j^ ^ reduced 
^rcd ~ and a gaussian scatter of the 19 GRBs around this 
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Fig. 2. Constraints on the cosmological parameters flu, i^A 
obtained with the updated sample of 19 GRBs presented in Tab. 
1 of N05 (to which GRB 051022 was added) in the wind den- 
sity profile case (WM). The solid (red) contours, obtained with 
the 19 GRBs alone, represent the 68.3%, 90% and 99% con- 
fidence regions on the pair of cosmological parameters. The 
contours obtained with the 156 SN la of the "Gold" sample of 
Riess et al. 2004 are shown by the long-dashed (blue) lines. 
The joint GRBh-SN constraints are represented by the shaded 
contours. We also report (dashed-red contours) the constraints 
obtained by assuming a fixed (i.e. cosmology invariant) hnear 
correlation E'pcak « E^. 

around this correlation with a a = 0.08) in the HM and WM 
case, respectively. 

3. Cosmological Constraints 

To the aim of constraining the cosmological parameters f^M 
and Ha we use the E-pcak-E^ correlations following the 
Bayesian method proposed by Firmani et al. (2005). The power 
of this optimized method is to circumvent the "circularity prob- 
lem" arising from the fact that these correlations cannot be cal- 
ibrated with the present sample of GRBs (see also Sec. 4). 

3. 1. Homogeneous density 

First we assume a HM and, through Eq. 1, derive the cosmo- 
logical constraints reported in Fig.^ The results, obtained with 
the present sample of 19 GRBs, are fully consistent with those 
obtained by Firmani et al. 2005 with the smaller (and slightly 
different) sample of 15 events (Tab. 1 of GGL04). The slightly 
larger contours reported in Fig. ^ are due to the changes of 
the parameters (and associated uncertainties) for GRB 991216, 
011211, 020124, 020405, 020813, 030226, 030328, 030329, 
030429 present in both samples (see N05 for a detailed discus- 



sion). The joint GRBh-SN fit (shaded regions in Fig.Q is dom- 
inated by the small contours determined with the 156 SN la 
of the "Gold" sample of Riess et al. (2004). However, we note 
that the GRBh-SN la joint fit is consistent at the 68% confidence 
level with the concordance model r^M, — (0.3, 0.7). 

3.2. Wind density profile 

In the case of a wind density profile we can still derive the 
constraints on the cosmological parameters adopting Eq. 2 to 
compute the jet opening angle. We make the simplest assump- 
tion of A* = 1. This choice corresponds to assuming a typ- 
ical mass loss rate A/^ — 10^^ Mq yr^^ and wind velocity 
Ww = 10'^ km s~^. We do not have any knowledge of the uncer- 
tainty associated to the parameter A*. However, as we want to 
compare these results with those obtained in the HM case (sec. 
3.1) where a typical uncertainty was asssociated to the den- 
sity parameter n, we assumed a 20% error on A^. The slightly 
smaller scatter of the Ghirlanda correlation in the case of the 
wind density profile and its better reduced — 1-14 (for 17 
dof), compared to the homogeneous case, are responsible for 
the more stringent constraints on the cosmological parameters 
in the ft plane. Our results are reported in Fig.|2] The 68% con- 
findence contours obtained with GRBs are still consistent with 
the concordance model although the center of the contours cor- 
responds to a quite large value of 17 a ■ 

On the other hand, in the WM case the Ep^ak-E^ correla- 
tion is linear (in the standard cosmology), and thus it has the 
remarkable property to be "Lorentz invariant, since both the 
E^ and E'pcak are boosted by ^ 2r when transforming from 
the comoving to the observer frame if we assume to view a 
uniform jet within the cone defined by its aperture angle. This 
property makes this correlation easier to be interpreted theoret- 
ically (see N05). We have then repeated our calculation assum- 
ing that the correlation remains linear in any cosmology. As 
discussed in N05, the linear Ghirlanda correlation obtained in 
the wind case is still consistent with the empirical correlation 
found by LZ05 and might have important implication for the 
physical interpretations of GRBs. If we then fix the slope to 1, 
the only free parameter remaining is the normalization. In this 
case we can find even better constraints (dotted red contours in 
Fig.EJ. 

4. GRB sample simulation 

In order to fully appreciate the potential use of GRBs for the 
cosmological investigation, we simulate a sample of bursts 
comparable in number to the "Gold" sample of 156 SN la. 
Similar simulations have been presented in the literature (Xu, 
Dai & Liang 2005; Liang & Zhang 2005); however different 
a-priori assumptions can be made on the properties of the sim- 
ulated sample and the results are clearly dependent on these 
assumptions. In particular, simulations based on the observed 
parameters (Xu et al. 2005) strongly depend on the selection 
effects on these quantities. 

We adopt here a method which makes use of the intrin- 
sic properties of GRBs as described by the Ghirlanda and the 
Amati correlations (Amati et al. 2002). We use the most up- 
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dated version of these correlations (N05) as found with the 
sample of 19 GRBs. With this sample the Amati correlation 
results: 



100 keV 



(0.361 ±0.02) 



7.6 X 1052 erg 



0.57±0.02 



with a reduced x^cd = 5-22 for 16 dof. 

The assumptions of our simulated GRB sample are: 



(3) 



- we assume that GRBs have an "isotropic energy" function 
which is described by a powerlaw N{Eiso) oc E^^^ for 
£^iso,min < £^iso < -Biso.max (e.g. Firmani et al. 2004). 
Further we assume that GRBs follow the cosmic star for- 
mation rate (SFR); 

- we use the Amati correlation (Eq.|3ll to derive the peak en- 
ergy Spoak; 

- we model the scatter of the simulated GRBs around the 
Amati correlation with a gaussian distribution with ct = 0.3 
(which corresponds to the present scatter of the 18 GRBs 
around their best fit correlation - Eq.|3ji; 

- we use the Ghirlanda correlation as found with the 18 
GRBs (either Eq. 1 or Eq. 2) to calculate E^; 

- we model the scatter of the simulated GRBs around the 
Ghirlanda correlation with a gaussian distribution with cr = 
0.1 (0.08 for the wind case - see N05); 

- we derive the jet opening angle 9j and the corresponding jet 
break time tj ; 

- we assume that the simulated GRB spectra are described 
by a Band model spectrum (Band et al. 1993) with typical 
low and high energy spectral photon indices a = — 1.0 and 
f3 = —2.5 and require that the simulated GRB fluence in 
the 2-400 keV energy band is above a typical instrumental 
detection threshold of ^ 10^^ erg/cm^. This corresponds 
roughly to the present threshold of Hete-II in the same en- 
ergy band. 

Following the procedure described above we built a sample 
of 150 GRBs with the relevant parameters: z, Ei^o, Ep, t j . The 
errors associated to these parameters are assumed to be cos- 
mology invariant and they are set to 10%, 20% and 20% for 
£■130, Ep, f j respectively. For our simulation we considered the 
Ghirlanda correlation in the WM case and adopted the SFR#2 
of Porciani & Madau (2001). We modeled the GRB intrinsic 
isotropic energy function with a powerlaw with S = —1.3 be- 
tween two limiting energies (10*^-10^^ erg). This particular 
choice of parameters is due to the requirement that the dis- 
tributions of the relevant quantities (shown in Fig. 3) of the 
simulated sample are consistent with the same distributions for 
the present sample of 19 GRBs. The sample is simulated in 
the standard cosmology (Hm = 0.3 and f^A = h = 0.7). We 
show the distribution of z,tb,Ep and Eiso for the 150 simu- 
lated bursts in comparison with the same distributions of the 
19 GRBs in Fig. 3. We also note that by choosing a steeper 
energy function we obtain a much larger number of XRF and 
XRR with respect to normal GRBs. 

The results obtained with the sample of 150 simulated 
GRBs is presented in Fig. |4] in this case the constraints are 
comparable with those obtained with SN la. The minimum of 
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Fig. 3. Distributions of redshift z, jet break time tb, observed 
peak energy E°]^l^ and isotropic equivalent energy Eiso for the 
150 simulated bursts (hatched histograms). Also show (solid 
filled histogram) are the distributions of the sample of 18 GRBs 
of N05 used to derive the ii^pcak — Eiso correlation. 



the GRB contours (cross in Fig. 4) corresponds to 17 j\/ = 0.27 
and = 0.72. By comparing the la contours of GRB alone 
from Fig. 4 to the same contours (solid line) of Fig. 2 (obtained 
with the 19 GRBs), we note that there is an improvement (of 
roughly a factor 10) with the sample of 150 bursts. 

It is evident the different orientation of the GRB contours 
(see Ghisellini et al. 2005) with respect to SN la due to the 
"topology" of the luminosity distance as a function of the ^Im- 
r^A parameters. Most GRBs of our simulated sample are at 
z ^ 2, and this explains the tilt of the contours obtained in the 
plane (Ghisellini et al. 2005). Clearly the simulated sample 
depends on the assumptions: in particular we have no knowl- 
edge of the burst intrinsic energy function N{E^^iso)- However, 
if we accept the hypothesis to model it with a simple powerlaw, 
we can change the slope and also include the effect of the red- 
shift evolution. We tested the dependence of these assumptions 
on the constraints reported in Fig. 4 and found that the major 
effect of assuming different S values and a (1 + z) evolution- 
ary factor is to change the redshift distribution of the simulated 
sample and therefore to change the orientation of the GRB con- 
tours in Fig. 4. The same happens if we adopt, for the same 
choice of parameters reported above, a different SFR. 

Further, we can also use the CMB priors. First we assume 
the cosmological constant model with the 2 CMB priors, i.e 
(i) iltot = 1 and (ii) fl^ = 0.14/ h^. In this case the only 
free parameter is h (or equivalently ilM).We obtain the best fit 
values of = 0.27 ±0.02 and = 0.73 ± 0.02. 
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Fig. 4. Wind density profile case. Constraints on the cosmo- 
logical parameters 57m, i^A obtained with the simulated sample 
of 150 GRBs. The soUd (red) contours, obtained with the 19 
GRBs alone, represent the 68.3%, 90% and 99% confidence 
regions on the pair of cosmological parameters. The contours 
obtained with the 156 SN la of the "Gold" sample of Riess et al. 
2004 are shown by the dashed (blue) Unes. The joint GRB+SN 
constraints are represented by the shaded contours. 
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Fig. 5. Constraints on the cosmological parameters wq, flu 
obtained with the 150 simulted GRBs (red contours) compared 
with the same contours obtained with the 156 SNIa of the 
"Gold" sample. A flat universe is assumed (iltot = !)■ 



Fig. 6. Constraints on the cosmological parameters wo, Wa ob- 
tained with the 150 simulted GRBs (red contours) compared 
with the same contours obtained with the 156 SNIa of the 
"Gold" sample. A flat universe is assumed (fitot = !)■ 

4.1. Dark Energy EOS 

One of the major promises of the cosmological use of GRBs 
is related to the possibility to study the nature of Dark Energy 
with such a class of "standard candles" extending out to very 
large redshifts. With the present sample of 19 GRBs we can 
explore the equation of state (EOS) of DE, which can be 
parametrized in different ways. Given the already considerably 
large dispersion of GRB redshifts (i.e. between 0.168 to 3.2 
for the 19 GRBs of our sample) we adopt the parametrization 
proposed by Linder & Huterer (2005) for the EOS of DE, i.e. 
P = w{z)p, where: 

, , WaZ 

w{z) =wo + —— (4) 
1 + z 

With this assumption the luminosity distance, as derived from 
the Friedmann equations, is 

c(l + z) 

+ (1 - r!M)(l + z)3+3-o+3». (^_3^^_^yi/2 

(5) 

which depends on the (wo,Wa) parameters. Note that Eq.|5]is 
derived with the prior of a flat Universe. 

First we can assume the CMB prior of a flat universe to- 
gether with the assumption of a non-evolving equation of state 
of the Dark Energy (i.e. Wa = 0). We show in Fig. 5 the con- 
tours obtained with the sample of 150 simulated GRBs and 
compare with the same constraints derived with the 156 SN la 
of the "Gold" sample. The constraints on Wa, wq are reported 
in Fig. 6, assuming ftyi = 0.3. 
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Fig. 7. Calibration of the E'pcak-^-^ correlation. For differ- 
ent samples of GRBs (6, 12, 18 - corresponding to the solid, 
dotted and dashed lines respectively) we show the maximum 
variation Ag of the slope of the correlation for any cosmol- 
ogy G (0, 1.5) as a function of the redshift dispersion of the 
GRBs dz. The dot dashed line represents the Umit of variation 
of 1% of the slope of the correlation. Data points have been 
shifted along the abscissa for graphycal purposes. 



5. The calibration of tlie spectrai energy 
correlations 

The cosmologiciil use of the £^peak = K ■ correlation suf- 
fers from the so called "circularity problem" (Ghirlanda et al. 
2004b, Ghisellini et al. 2005): this means that both the slope g 
and the normalization K of the correlation are cosmology de- 
pendent. In fact, of the two rest frame quantities iJpeak and 
that are used to compute the Ghirlanda correlation, the second 
one (E^^ depends on the cosmological model through the lumi- 
nosity distance {z;^. In principle this issue could be solved 
(a) with a considerably large number of calibrators, i.e. low red- 
shift GRBs for which the luminosity distance dz, is practically 
independent from the cosmological parameters, or (b) with a 
convincing theoretical interpretation of the physical nature of 



this correlation. In both cases the slope of the correlation would 
be fixed. 

Case (a) could be realized with 5-6 GRBs at ^ < 0.1. 
However, if (long) GRBs are produced by the core-collapse of 
massive stars, their rate is mainly regulated by the cosmic SFR 
and, therefore, the probability of detecting events at z < 0.1 
is small (~ 2 x 10~^). This number should be convolved with 
the GRB luminosity function: with the assumptions (described 
in Sec. 4) of our simulation we estimate that ~ 1.3% of the 
150 GRBs should be at low redshifts (i.e. z < 0.4). Instead 
we should expect to have more chances to detect a consid- 
erable number (up to ~ 31%) of intermediate redshift GRBs 
(z ~ 1 — 2) where the cosmic SFR peaks. 

For this reason we explore the possibility to cahbrate the 
correlation using a sufficient number of GRBs within a small 
redshift bin centered around any redshift. In fact, if we could 
have a sample of GRBs all at the same redshift the slope of 
the Ghirlanda correlation would be cosmology independent. 
Our objective is, therefore, to estimate the minimum number 
of GRBs N within a redshift bin dz centered around a certain 
redshift Zc which are required to cahbrate the correlation. 

In practice the method consists in fitting the correlation for 

every choice of J7 using a set of N GRBs distributed in the 
interval dz (centered around Zc). We consider the correlation to 
be calibrated (i.e. its slope to be cosmology independent) if the 
change of the slope g is less than 1%. 

The free parameters of this test are the number of GRBs A'^, 
the "redshift sUce" dz and the central value of the redshift dis- 
tribution Zc- By Monte Carlo technique we use the same sam- 
ple simulated in Sec. 4 under the WM assumption to minimize 
the variation of Ag{fl; N, dz, Zc) over the VLm, & (0, 1.5) 
plane as a function of the free parameters {N,dz,Zc). 

We tested different values of Zc and different redshift dis- 
persions dz e (0.05, 0.5). We required a minimum number of 
6 GRBs to fit the correlation in order to have at least 4 degrees 
of freedom. We report our results in Fig. 7. We show the varia- 
tion of A (7 as a function of dz for different samples (6, 12, 18 
GRBs - solid, dotted and dashed curves in Fig. 7). The error 
bars show the width of the distribution of the simulation results 
and not the uncertainty on the average value. At any redshift 
the fewer the number A'' of GRBs the larger the change of Ag 
(for the same dz) because the correlation is less constrained. 
The dependence from Zc is instead different: for larger Zc we 
require a smaller bin dz to keep Ag small. 

From the curves reported in Fig. 7 we can conclude that al- 
ready 12 GRBs with z G (0.9, 1.1) might be used to cahbrate 
the slope of the Epcak-E-y correlation. At redshift Zc = 2 in- 
stead we require a smaller redshift bin i.e. z G (1.95, 1.05). 
We find that = 12 GRBs with z e (0.45, 0.75) can be used 
to achieve the same 1% precision in the calibration. However, 
one key ingredient is that the GRBs used to calibrate the corre- 
lation do not have the same peak energy otherwise they would 
collapse in one point in the Epcak-E^ plane. Within the present 
sample of 19 GRBs there are only 4 GRBs within the redshift 
interval 0.4-0.8 (i.e. 050525, 041006, 020405 and 051022) and 
2 of these (050525 and 041006) have a very similar i^peak- 
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6. Conclusions 

We have presented the cosmological constraints that can be ob- 
tained with the present sample of 19 GRBs for which all the 
relevant quantities, i.e. redshift z, peak energy -Bpcak and jet 
break time t\y, has been properly estimated and published in the 
literature. 

Following the results of N05, where the E^q^\^-E^ corre- 
lation found under the hypothesis of a uniform jet model with 
either a homogeneous (HM) or a wind density profile (WM) 
were presented, we derived the constraints on 51m and 57a in 
these two scenarios. Tighter constraints are obtained in the WM 
model. 

We also presented the future of GRBs as cosmological 
tools. By simulating a sample of 150 GRBs (which is com- 
parable in number to the present sample of "Gold" SN la) we 
showed that similar (to SN la) tight constraints can be obtained 
with GRBs either on the present universe composition and on 
the nature and evolution of DE (parametrized with and Wa)- 
We remark that the collection of such a large sample of GRBs 
which can be used as standard candles requires an accurate 
measurement of their prompt and afterglow properties. In par- 
ticular a wide energy spectral coverage is required to constrain 
the peak energy and properly compute the bolometric corrected 
GRB energetics. 

Finally, a large sample of GRBs would help in calibrating 
the slope of the E^c&v-E^ correlation (in both the HM or WM 
case). To this aim one would require to collect enough low 
redshift GRBs. However, we can relax this requirement: we 
showed that even if not concentrated at very low z, it is suf- 
ficient to have a dozen GRBs with a similar redshift to find the 
slope of the correlation in a cosmology-independent way at the 
level of 1 % accuracy. 
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